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Abstract 



The multichannel Kondo model with SU(A^) spin symmetry and SU(i^) chan- 
nel symmetry is considered. The impurity spin is chosen to transform as an 
antisymmetric representation of SU(A^), corresponding to a fixed number of 
Abrikosov fermions J2afafa = Q- more than one channel {K > 1), and 
all values of N and Q, the model displays non- Fermi behaviour associated 
with the overscreening of the impurity spin. Universal low-temperature ther- 
modynamic and transport properties of this non-Fermi liquid state are com- 
puted using conformal field theory methods. A large- limit of the model is 
then considered, in which K/N = 7 and Q/N = qq are held fixed. Spectral 
densities satisfy coupled integral equations in this limit, corresponding to a 
(time-dependent) saddle-point. A low frequency, low-temperature analysis of 
these equations reveals universal scaling properties in the variable co/T, in 
agreement with conformal invariance. The universal scaling form is obtained 
analytically and used to compute the low-temperature universal properties of 
the model in the large-N limit, such as the T = residual entropy and resid- 
ual resistivity, and the critical exponents associated with the specific heat and 
susceptibility. The connections with the "non-crossing approximation" and 
the previous work of Cox and Ruckenstein are discussed. 

PACS numbers: 75.20 Hr, 75.30 Mb, 71.10 Hf. Preprint LPTENS : 97/55 
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I. INTRODUCTION AND MODEL 



Multichannel Kondo impurity models have recently attracted considerable attention, 
for several reasons. First, in the overscreened case, they provide an explicit example of a non- 
Fermi liquid ground-state. Second, these models can be studied by a variety of controlled 
techniques, and provide invaluable testing grounds for theoretical methods dealing with 
correlated electron systems. One of the most recent and fruitful development in this respect 
has been the conformal field-theory approach developed by Affleck and Ludwig . Finally, 
multichannel models have experimental relevance to tunneling phenomena in quantum dots 
and two level systems 0, and possibly also to some heavy-fermion compounds 0. 
In this paper, we consider a generalisation of the multi-channel Kondo model, in which the 
spin symmetry group is extended from SU(2) to SU(A^). In addition, the model has a SU(ii') 
symmetry among the K "channels" (flavours) of conduction electrons. We focus here on 
spin representations which are such that the model is in the non-Fermi liquid overscreened 
regime. We shall derive in this paper several universal properties of this SU(A^) ® SU(i^) 
Kondo model in the low-temperature regime. Speciflcally, we obtain: the zero-temperature 
residual entropy, the zero-temperature impurity resistivity and T-matrix, and the critical 
exponents governing the leading low-temperature behaviour of the impurity speciflc heat, 
susceptibility and resistivity. 

These results will be obtained using two different approaches. It is one of the main moti- 
vations of this paper to compare these two approaches in some detail. First (Sec.pl), we 
apply CFT methods to study the model for general values of and K. Then, we study 
the limit of large and K, with K/N = 7 flxed. This limit was previously considered 
by Cox and Ruckenstein 0, in connection with the "non-crossing approximation" (NCA) . 
There is a crucial difference between our work and that of Ref. however, which is that 
we keep track of the quantum number specifying the spin representation of the impurity 
by imposing a constraint (on the Abrikosov fermions representing the impurity) which also 
scales proportionally to N. As a result, the solution of the model at large N follows from a 
true saddle-point principle, with controllable fluctuations in Hence, a detailed quan- 

titative comparison of the large-N limit to the CFT results can be made. The saddle-point 
equations are coupled integral equations similar in structure to those of the NCA, except for 
the different handling of the constraint. The T = impurity entropy and residual resistivity 
are obtained for the flrst time in analytical form in this paper from a low-energy analysis 
of these coupled integral equations and shown to agree with the large-N limit of the CFT 
results. We also demonstrate that the spectral functions resulting from these equations take 
a universal scaling form in the limit uj,T 0, which is precisely that expected from the 
conformal invariance of the problem. 

The Hamiltonian of the model considered in this paper reads: 

K N N^-l 
p> i=l a=l A=l ^p'ia/3 

In this expression, ct^^ creates an electron in the conduction band, with momentum p, chan- 
nel (flavour) index i = 1, - ■ ■ ,K and SU(A^) spin index a = 1, - ■ ■ ,N. The conduction 
electrons transform under the fundamental representation of the SU(A^) group, with gener- 
ators t^^ {A = 1, ■ ■ ■ , N"^ — 1). They interact with a localised spin degree of freedom placed 
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at the origin, S = {S^, A = 1, ■ ■ ■ , N'^ — 1} which is assumed to transform under a given 
irreducible representation R of the SU(iV) group. 

In the one-channel case [K = 1), and when R is taken to be the fundamental representation, 
this is the Coqblin-Schrieffer model of a conduction gas interacting with a localised atomic 
level with angular momentum j {N = 2j ' + 1) [^. In this article, we are interested in 
the possible non-Fermi liquid behaviour associated with the multi-channel generalisation 
{K > 1) [Q 0. We shall mostly focus on the case where R corresponds to antisymmetric 
tensors of Q indices, i.e the Young tableau associated with R is made of a single column 
of Q indices. In that case, it is convenient to use an explicit representation of the localised 
spin in terms of N species of auxiliary fermions fa (a = 1, ■ ■ ■ , N), constrained to obey: 



N 



E /I/- = Q (2) 
«=i 

so that the N'^ — 1 (traceless) components of S can be represented as: Saf3 = fafp ~ ^^a(3- 
For these choices of i?, the Hamiltonian can be written as (after a reshuffling of indices using 
a Fierz identity ): 



K N 



p i=l Q=l jjp'iap 

In a recent paper |TD|, the case of a symmetric representation of the impurity spin (corre- 



sponding to a Young tableau made of a single line of P boxes) has been considered by two 
of us. In that case, a transition from overscreening to underscreening is found as a func- 
tion of the "size" P of the impurity spin. In contrast, the antisymmetric representations 
considered in the present paper always lead to overscreening (except for K = 1 which is 
exactly screened), as shown below. As long as only the overscreened regime is considered, 
the analysis of the present paper applies to symmetric representations as well, up to some 
straightforward replacements. 



II. STRONG COUPLING ANALYSIS 

It is easily checked that a weak antiferromagnetic coupling {Jk > 0) grows under renormal- 
isation for all K and A^, and all representations R of the local spin. What is needed is a 
physical argument in order to determine whether the R.G flow takes Jk all the way to strong 
coupling (underscreened or exactly screened cases), or whether an intermediate non- Fermi 
liquid flxed point exists (overscreened cases). 

Following the Nozieres and Blandin analysis of the multichannel SU(2) model, we con- 
sider the strong-coupling flxed point Jk = +00. In this limit, the impurity spin binds a 
certain number of conduction electrons (at most NK because of the Pauli principle). The 
resulting bound-state corresponds to a new spin representation Rgc which is dictated by the 
minimisation of the Kondo energy. For the speciflc choices of R above, we have proven that: 

• In the one-channel case {K = 1) and for arbitrary A^ and Q, Rsc is the "singlet" 
representation (of dimension d{Rsc) = !)• It is obtained by binding N — Q conduction 
electrons to the Q pseudo-fermions /„. The impurity spin is thus exactly screened at 
strong-coupling. 
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• For all mult i- channel cases {K > 2, arbitrary and Q), the ground-state at the 
strong-coupling fixed point is the representation Rgc characterised by a rectangular 
Young tableau with N — Q lines and K — 1 columns. Its dimension d{Rsc) {i-e the 
degeneracy of the strong-coupling bound-state) is larger than the degeneracy at zero 
coupling, given by d{R) = (J) = N\/Q\{N -Q)\. 

The Young-tableau associated with the strong-coupling state in both cases is depicted in 
Fig.|l]. The detailed proof of these statements and the explicit construction of Rgc are given 
in Appendix These properties are sufficient to establish that : 

• In the one-channel case, the strong coupling fixed point is stable under R.G., and 
hence the impurity spin is exactly screened by the Kondo effect. 

• In the multi-channel direct R.G flow from weak to strong-coupling is impossible, 
thereby suggesting the existence of an intermediate coupling flxed-point ( "overscreen- 
ing"). 

The connection between these statements and the above results on the nature and degeneracy 
of the strong-coupling bound-state is clear on physical grounds. Indeed, it is not possible 
to flow under renormalisation from a flxed point with a lower ground-state degeneracy to a 
flxed-point with a higher one, because the effective number of degrees of freedom can only 
decrease under R.G. Hence, no flow away from the strong-coupling flxed point is possible 
in the one-channel {K = 1) case since the strong-coupling state is nondegenerate. Also, no 
direct flow from weak to strong coupling is possible foYK>2 since d{Rsc) > d{R). These 
statements can be made more rigorous by considering the impurity entropy deflned as: 

Simp = lim lim [S{T) - SbuikiT)] (4) 

T — ^0 V — ^oo 

Where Sbuik denotes the contribution to the entropy which is proportional to the volume 
V (and is simply the contribution of the conduction electron gas), and care has been taken 
in specifying the order of the inflnite-volume and zero-temperature limits. At the weak 
coupling fixed point Simp{JK = 0) = lnd{R), while Simp{JK = co) = hid{Rsc) at strong- 
coupling. Simp niust decrease under renormalisation 0, a property which is the analogue 
for boundary critical phenomena to Zamolodchikov's "c-theorem" in the bulk. This suggests 
a R.G fiow of the kind indicated above. This conclusion can of course be confirmed by a 
perturbative calculation (in the hopping amplitude) around the strong-couphng fixed-point 

i- 

The value of Simp will be calculated below at the intermediate fixed-point in the overscreened 
case, and found to be non-integer (as in the N = 2 case 0). 

III. CONFORMAL FIELD THEORY APPROACH 

Having established the existence of an intermediate fixed-point for K > 2, we sketch some 
of its properties that can be obtained from conformal field-theory (CFT) methods. This is 
a straightforward extension to the SU(A^) case of Affleck and Ludwig's approach for SU(2) 
IQJ^. The aim of this section is not to present a complete conformal field-theory solution, but 
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simply to derive those properties which will be compared with the large-N explicit solution 
given below. 

In the CFT approach, the model (|1]) is first mapped at low-energy onto a 1 + 1-dimensional 
model of NK chiral fermions. At a fixed point, this model has a local conformal symmetry 
based on the Kac-Moody algebra SUk{N)s C?> SU N{K)f (g) f/(l)c corresponding to the spin- 
flavour- charge decomposition of the degrees of freedom. The free-fermion spectrum at the 
weak-coupling fixed point can be organised in multiplets of this symmetry algebra: to each 
level corresponds a primary operator in the spin, flavour and charge sectors. A major insight 
1^ is then that the spectrum at the infra-red stable, intermediate coupling fixed point can 
be obtained from a "fusion principle". Specifically, the spectrum is obtained by acting, in 
the spin sector, on the primary operator associated with a given free-fermion state, with 
the primary operator of the SUk{N)s algebra corresponding to the representation R of the 
impurity spin (leaving unchanged the flavour and charge sectors). The "fusion rules" of the 
algebra determine the new operators associated with each energy level at the intermediate 
fixed point. 

This fusion principle also relates the impurity entropy Simp as defined above to the "modular 
S- matrix" Sq of the SUk{N) algebra (this is the matrix which specifies the action of a 
modular transformation on the irreducible characters of the algebra corresponding to a 
given irreducible representation R). Specifically, denoting by i? = the trivial (identity) 
representation: 

Simp — ^^^Q (^) 



The expression of the modular S-matrix for SUk{N) can be found in the litterature |]TT|. We 



have found particularly useful to make use of an elegant formulation introduced by Douglas 
I^ , which is briefly explained in Appendices ^ and |^. For the representations R associated 



to a single column of length Q (corresponding to (||)), one finds using this representation: 

fi- = n —zit- (6) 

n=l N+K 

It is easily checked that indeed Simp < ln(i(i?) for all values of A^, Q and K. Note also that 
this expression correctly yields Simp = in the exactly screened case K = 1 (for arbitrary 
N,Q). 

The low-temperature behaviour of various physical quantities can also be obtained from the 
CFT approach. At the intermediate coupling fixed point, the local impurity spin acquires 
the scaling dimension of the primary operator of the SUk{N)s algebra associated with the 
(A^^ — 1 dimensional) adjoint representation of SU(A^). Its conformal dimension (such 
that < S{0)S{t) >~ at T = 0) reads: 

^' - nTk pi 

Integrating this correlation function, this implies that the local susceptibility xioc oc J^^j'^ < 
S{0)S{t) > dr (corresponding to the coupling of an external field to the impurity spin only) 
diverges at low temperature when K > N, while it remains finite for K < N: 
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K >N : xioc ~ (^j 

K = N : x;oc~lnl/T 

K < N : xioc ~ const. (8) 

Exactly at the fixed point, the singular contributions to the specific heat and impurity sus- 
ceptibility Ximp = X ~ Xbulk (defined by coupling a magnetic field to the total spin density) 
vanish |3,[T3| . Indeed, the singular behaviour is controlled by the leading irrelevant operator 
compatible with all symmetries that can be generated. An obvious candidate for this opera- 
tor is the spin, fiavour and charge singlet obtained by contracting the spin current with the 
adjoint primary operator above, J_i ■ 0. It has dimension 1 + = 1 + A^. This leads to 
a singular contribution to the impurity susceptibility (arising from perturbation theory at 
second order in the irrelevant operator P,p!3|) of the same nature than for xioc given in (^: 

Ximp ~ Xioc- 

Another irrelevant operator can be constructed in the flavour sector in an analogous manner, 
namely jLi ■ (pf . This operator has dimension K/{N + K)^ which is thus lower than the 
above operator in the spin sector when K < N . This operator could a priori contribute to 
the low-temperature behaviour of the specific heat, which would lead to a divergent specific 
heat coefficient C/T for both K > N and K < N. On the basis of the explicit large- 
N calculation presented below, we believe however that this operator is not generated for 
model (0) when the conduction band d.o.s is taken to be perfectly fiat and the cutoff is 
taken to infinity (conformal limit), so that the specific heat ratio has the same behaviour 
than the susceptibilities above: C/T ~ xioc ~ Ximp- If the model is extended to an impurity 
spin with internal fiavour degrees of freedom, this operator will however show up (leading 
to C/T ~ j'~{N-K)/{N+K) ^ ^ jY^_ g^jgQ appgarg (^ggg below) if an Anderson model 

generalisation of model (|I]) is considered away from particle-hole symmetry. 



IV. SADDLE-POINT EQUATIONS IN THE LARGE-N LIMIT 

We now turn to the analysis of the large-N limit of this model. This will be done by setting: 

K = , Jk =^ (9) 

and taking the limit N ^ oo for fixed values of 7, and J, so that the number of channels is 
also taken to be large. In Ref. |^ (see also 0), Cox and Ruckenstein considered this limit 
while holding Q fixed [Q = 1). They obtained in this limit identical results to those of the 
"non-crossing approximation" (NCA) [^]. Here, we shall proceed in a different manner by 
taking Q to be large as well: 

Q = qoN (10) 

This insures that a true saddle-point exists, with controllable fluctuations order by order in 
1/A^. It will also allow us to study the dependence on the representation R of the local spin, 
parametrised by go The approach of Ref. [0 is recovered in the limit go (or 1). 



The action corresponding to the functional integral formulation of model (^) reads: 
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+ !^ dr Ea fi{r)drfair) + !^ dr ifi{T) (Ea Pair) fair) - qoN) 

+ i !^ dr Y^iap cl{r)cip{T) {fl{T)UT) - qoSap) (11) 

In this expression, the conduction electrons have been integrated out in the bulk, keeping 
only degrees of freedom at the impurity site. Go(^'^n) = Epl/l^i-^n — ^p) is the on-site 
Green's function associated with the conduction electron bath. In order to decouple the 
Kondo interaction, an auxiliary bosonic field Bi^r) is introduced in each channel, and the 
conduction electrons can be integrated out, leaving us with the effective action: 

Seff = /o^rfTEa/l(r)9,/«(r)+/o^dri/i(r)(Ea/l(T)/a(T)-goA^) + yo^^iTE.SlS« 

+ J! !i dr dr' Eia B,{T)fl{T)Go{T - t')B\{t')Ut') (12) 

The quartic term in this expression can be decoupled formally using two bi-local fields 
Q{t,t') and Q{t,t') conjugate to Y,iB\{T')Bi{T) and J2a fa{^)fa{T') respectively, leading 
to the action : 

f dTY.fi{r)drUT) + \ f drY^BlB.^ f dTiix{T)(£fl{T)UT)-qoN) 
Jo ^ J Jo Y ■'^ a 

-N JJdTdT'Q{r,r')G^\r-T')Q{r,r')- jj drdr' Q{r',T)^B+{r)B,{T') 

i 

- [[ dTdT'QiT,T')j:f^iT)UT') (13) 

The B and / fields can now be integrated out to yield: 

-N JJ dTdr'Q{r,T')G^\r-T')Q{r,r')-NqoJ drii^ir) 
-NTr ln[{-dr - ^/^(r))5(T - r') + g(r, r')] + i^Tr ln[^5(r - r') - Q{t\ r)] (14) 

This final form of the action involves only the three fields Q, Q and and scales globally as 
thanks to the scalings K = 'fN and Q = qoN. Hence, it can be solved by the saddle-point 
method in the large-N limit. At the saddle-point, /isp(T) = i\ becomes static and purely 
imaginary, while Qgp = Q{t — r') and Q^^ = Q{t — t') retain time- dependence but depend 
only on the time difference t — r' (they identify with the bosonic and fermionic self-energies, 
respectively) . 

The final form of the coupled saddle-point equations for the fermionic and bosonic Green's 
functions G/(r) = - < Tf{T)p{0) >, GbIt) =< TB{t)B\Q) > and for the Lagrange 
multiplier field read: 

E^(r) = ^G,{t)Gs{t) , EB(r) = Go(r)G/(r) (15) 
where the self-energies and are defined by: 

Gj\iuJn) =iuJn + >^- ^f{i(^n) , G^'^ (iUn) = j - T,B{il^n) (16) 
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In these expressions Un = {^n + and z/„ = 2mi/P denote fermionic and bosonic 

Matsubara frequencies. Let us note that the field B{t) is simply a commuting auxiliary 
field rather than a true boson (its equal-time commutator vanishes). As a result Gb{t) is 
a /3— periodic function, but does not share the usual other properties of a bosonic Green 
function (in particular at high frequency) . 
Finally, A is determined by the third saddle point equation : 

Gf{r = 0-) = i 5: G/(z^„)e^""°" = go (17) 

P n 

These equations are identical in structure to the usual NCA equations, except for the last 



equation ([T7|) which implements the constraint and allows us to keep track of the choice of 



representation for the impurity spin. 



V. SCALING ANALYSIS AT LOW FREQUENCY AND TEMPERATURE 

A. General considerations 



The analysis of the NCA equations in Ref. |T5[ can be applied in order to find the behaviour 
of the Green's functions in the low temperature, long time regime defined by: T^^ <^ t <^ 
(3 —* oo (where Tr- is the Kondo temperature). In this regime, a power-law decay of the 
Green's functions is found: 

Gfir) ~ , Gnir) ~ ^ , (T^^ « r « /3 - oo) (18) 

The scaling dimensions 2Aj and 2Ab can be determined explicitly by inserting this form 
into the above saddle-point equations and making a low-frequency analysis, as explained in 
Appendix y. This yields: 

2A/ = , 2Ab = (19) 

1+7 1+7 

The overall consistency of (|1^,0) at large time also constrains the product of amplitudes 
AjAb (Eq. (|C^ ) in Appendix 0) and dictates the behaviour of the self-energies (denoting by 
po = — ImG'o(^0^)/7r the conduction bath density of states at the Fermi level) : 

S^r) ~7Po^B (-j , SB(r) ~ A;po (-j (20) 



together with the sum rule: 



SB(cc^ = 0,/5 = oo) = i (21) 



The expression (|T9|) of the scaling dimensions and is in complete agreement with 



the CFT result. Indeed, the fermionic field transforms as the fundamental representation of 
the SU{N)k spin algebra, while the auxiliary bosonic field transforms as the fundamental 
representation of the SU{K)j^ fiavour algebra, leading to: 2Aj = ^^^"^^ and 2Ab = 
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K{N+K) ' which agree with (|T^) in the large-N hmit. Also, the local impurity spin correlation 
function, given in the large-N limit by (S'(0)5'(r)) oc G/(r)Gj(— r) is found to decay as 

^2As „,UT. A _ 2A 



with = 2A/ = 1/(1 + 7) in agreement with the CFT result (0). 
As will be shown below however, these asymptotic behaviours at T = do not provide 
enough information to allow for the computation of the low-temperature behaviour of the 
impurity free-energy and to determine the T = impurity entropy (H). One actually needs 
to determine the Green's functions in the limit t, f3 ^ 00, but for an arbitrary value of the 
ratio t/P {i.e to analyze the low-temperature, low-frequency behaviour of (15,16) keeping 
the ratio u/T fixed [|T7|). It is easily seen that in this limit the Green's functions and their 
associated spectral densities pj b(u;) = —^ImGf^Bi^ + ^0^) obey a scaling behaviour : 



-2A 



f.B 



9f,B 




<ti T , P ; t/P arbitrary) 



(22a) 



Pfiu;) = AfT'^f-' 4>f (^) , p^(^) = AsT'^--' <Pb (^) (22b) 

In these expressions, gf^B and are universal scaling functions which depend only on 7 
and go and not on the specific shape of the conduction hand or the cutoff. These scaling 
functions will now be found in explicit form. 



B. The particle-hole symmetric representation go = 1/2 

We shall first discuss the case where the representation R has Q = N/2 boxes (go = 1/2), 
for which there is a particle- hole symmetry among pseudo-fermions under f^ ^ f^. The 
expression of the scaling functions gf^B in that case can be easily guessed from general 
principles of conformal invariance. The idea is that, in the limit T^^ -C P with t/P fixed, 
the finite-temperature Green's function can be obtained from the T = Green's function by 
applying the conformal transformation z = exp{i2'n'T / P) Applying this to the T = 
power-law decay given by (0), one obtains the well-known result for the scaling functions 

{f^r/py. 



^7;(f;go = l/2) = - ^ ^ , gB{f;qo = l/2) = -(-. ^ (23) 

VsmTrr/ VsmTrr/ 

with the periodicity requirements gf{f + 1) = —gf{f), gBij + 1) = gBij)- Note that these 
functions satisfy the additional symmetry gf B{^ — t) = gf^B{^) indicating that they can 
only apply to the particle-hole symmetric case go = 1/2. The corresponding form of the 
scaling functions associated with the spectral densities ( |22b[ ) reads, with Cj = uj/T (after a 
calculation detailed in Appendix 0): 

0;(^, go = 1/2) = i(2vr)2A/-i cosh f "^^^^V|a,^^"^^^ 

M^. qo = 1/2) = i(27r)2A--i sinh f "^^"%f I^^""^^^ (24) 
Note that the u; — > singularity of the T = case is now recovered for uj ^ T: 
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1 / 1 \ 1-2A/.B 

go = 1/2) , ~ . ^ (25) 

It is an interesting calculation, performed in detail in Appendix y, to check that indeed these 
scaling functions do solve the large- equations ([15| , [T6|) at finite temperature in the scaling 
regime. That NCA-like integral equations obey the finite-temperature scaling properties 
dictated by conformal invariance has not, to our knowledge, been pointed out in the previous 
litterature. 



C. General values of go 

1. Spectral asymmetry 

Let us move to the general case of representations with ^ in which the particle-hole 
symmetry between pseudo-fermions is broken. The exponent of the power-law singularity 
in the T = spectral densities is not affected by this asymmetry. It does induce however 
an asymmetry of the prefactors associated with positive and negative frequencies as ^ 0. 
We introduce an angle 6 to parametrize this asymmetry, defined such that: 

a-N ,/ sinfTrAf + 0) , x sinfvrAf — 

pj{u 0+) ~ M7, 9) IJ^^. ^ , P/(^ - 0-) ~ M7, 0) (_^)/-2A, (26) 

where h{'y, 0) is a constant prefactor. The explicit dependence of on go will be derived 
below. This corresponds to the following analytic behaviour of the Green's function in the 
complex frequency plane, as 2; — > 0: 

Gf{z)r^h{j,9) lmz>0 (27) 

Equivalently, this means that the symmetry Gf{(3 — t) = Gf{T) is broken, and that the 
scaling function Qfir) must satisfy (from the behaviour of its Fourier transform): 

^ sin(7r A; + 9) 
gf{l~) sm{nAf-e) ^ ^ 

We have found, by an explicit analysis of the saddle-point and constraint equations in the 
scaling regime, which is detailed in Appendix 0, that the full scaling functions for this 
asymmetric case are very simply related to the symmetric ones at go = 1/2, through: 

where the parameter a is simply related to 6 so as to obey (pS]) : 

« = ln— )^il±2^ J (30) 
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Fourier transforming, this leads to the scahng functions for the spectral densities: 



cosh - 

= cosh ^ cosh f '^^^^ + = 
sinh ^ 



^b{^, go) = ^ (t>B{^ + a,qo = 1/2) (31) 



sinh cosh f 



The thermal scaling function for the fermionic spectral density (pf and the bosonic one (pb 
are plotted in Fig.(|]) for various values of the asymmetry parameter a. We also note the 
expression for the maximally asymmetric case a —>■ —oo (corresponding, as shown below, to 
go — * 0, i.e. to the limit Q ^ as in the usual NCA). 



, . (2vr)^^^-^r(A, + .f)r(A,-^f) 

(hf(uj-a) ~ 62 , , ( — ^ ^ (32) 

^^^-^ 7rr(2A/) ^ ^ 



a— >— oo 



We also note, for further use, the expressions of the full Green's functions in the complex 
frequency plane (defined by gj/siz, a) = dcD^^^j^i^), in the scaling regime for Im^; > 0: 



gf{z,a) 



2z(27r) 



cosh fr(2A/) sin 27rA/ 

2n J \ ^ 27r 



z + a\ ^ / ^ z + a\ ( ^ ia\ z + a\ ,, 

A/ + ^^^ r A; -z— — cos ttA; - sm ttA; + z^— (33) 



2(27r)2^s-i 
gB{z,a) - -coghfr(2AB)sin27rAB 

r [^B + ^^^) r (^Ab - j sm (^ttAb - y j sm (^ttAb + j (34) 

The reader interested in the details of these calculations is directed to Appendix |C[ 
At this stage, the point which remains to be clarified is the explicit relation between the 
asymmetry parameter Q and the parameter go specifying the representation. This is the 
subject of the next section. 

Before turning to this point, we briefly comment on the CFT interpretation of the asymmetry 
parameter Q (or a) associated with the particle-hole asymmetry of the fermionic fields. The 
form ( [23| ) of the correlation functions at finite temperature in the scaling limit can be viewed 
as those of the exponential of a compact bosonic field with periodic boundary conditions. 
The asymmetric generalization ( p9D corresponds to a shifted boundary condition on the 
boson [i.e to a twisted boundary condition for its exponential). 



2. Relation between qq and 9 

Let us clarify the relation between the spectral asymmetry parameter 6, and the parame- 
ter go specifying the spin representation. That such a relation exists in universal form is 
a remarkable fact: indeed 6' is a low-energy parameter associated with the low-frequency 
behaviour of the spectral density, while go is the total pseudo-fermion number related by 
the constraint equation (|l^) to an integral of the spectral density over all frequencies. The 
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situation is similar to that of the Friedel sum rule in impurity models, or to Luttinger the- 
orem in a Fermi liquid, and indeed the derivation of the relation between go and 9 follows 
similar lines [^]. It is in a sense a Friedel sum rule for the quasiparticles carrying the spin 
degrees of freedom (namely, the pseudofermions fa)- 

We start from the constraint equation (O) written at zero-temperature as: 



— i lim ( — Gr{uj)e 



iijjt 



(35) 



In this expression, and below in this section, Gf{uj) and Gb{oj) denote the (Feynman) T = 
Green's functions while the retarded Green's functions are denoted by G^. Using analj^ic 
continuation of (0), we have : 



^/H 

(9 In GBi^ 



duj 

so that (|35|) can be rewritten as: 



duo 
duj 



-GAu) 



(36) 
(37) 



go 



i lim / — 

t^o+ J 27r 



dlxiG 



.00 



duo 



G, 



9Sy 



[uo 



[00 



duo 



7 



'd\nGB{uj) 
duo 



Gb{uj 



dT.B{uo) 
duo 



(38) 



In this expression, the bosonic part (which vanishes altogether) has been included in order to 
transform further the terms involving derivatives of the self energy, using analyticity. This 
transformation is only possible if both fermionic and bosonic terms are considered. This is 



because the Luttinger- Ward functional |T9[ of this model involves both Green's functions. 
It has a simple explicit expression which reads: 



^Lw{Gf,a,GB,) = E/ dtGo{t)Gf^a{-t)GBAt) 



(39) 



such that the saddle-point equations (|15|) are recovered by derivation: 



6G 



6G 



B,: 



-t) 



(40) 



From the existence of ^lw, we obtain the sum rule: 



/ dto S/(ci;) 

J — CXI \ 







(41) 



(Note that there is no logarithmic divergence in this integral). After integrating by parts 
and using the asymptotic behaviour of the two Green's functions in order to eliminate the 
boundary terms, we get : 



roo du f dlnGfjio) _ d\nGBiuj) \ 
^° ^ ^oo 27r I duo ^ duo j ^ 



(42) 
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Since G{uj) = G^{uj) for > and G{uj) = G^{uj) for < (with G^ the retarded Green's 
function), this can be transformed using:. 



- dudlnGfJu) , r dudlnGfJcu) , po du d f Gf^ 



27r duj J-oQ 2n dui J-oo 27i dui \ Gf^^[uj 

(43) 

The first integrals in the right hand side can be deformed in the upper plane and their sum 
vanishes |jT9|. Thus we obtain (denoting by argG the argument of G) : 

vrgo = argGf (0-) - argCf (-oo) - 7 (argCf (0") - argCK-oo)) (44) 

argGj(0~) directly follows from the parametrisation ( pTf ) defining 9. It can also be read off 
from the behaviour of the scaling function gf{z) for z = ±00. Thus, from (|33|) we can also 
read off arg 6*1^(0"): 

argGf (0") = ttAj - - TT , argGf (0") = 6 - nAf (45) 

Taking into account that Gf^cu) ~a;^-oo V"^ that ImG^ < we have argG^ (— 00) = 
— vr. Similarly, we have argG|^(— 00) = 0. Inserting these expressions into (^^, we finally 
obtain the desired relation between go and 9 (or a): 

1+7 

This, together with (|33[) , fully determines the universal scaling form of the spectral functions 
in the low- frequency, low temperature limit. 



VI. PHYSICAL QUANTITIES AND COMPARISON WITH THE CFT 

APPROACH 

A. Impurity residual entropy at T = 

The impurity contribution to the free-energy (per colour of spin) fi^p = {F — Ff,uik) / N reads, 
at the saddle point: 

fimp = qoX + Tj2\nGf{tuj^^) - ^Tj2\nGB{tiyn) - f drT.f{r)Gf{-T) (47) 

This expression can be derived either directly from the saddle-point effective action (in 
which case the last term arises from the quadratic term in Q and Q), or from the relation 
between the free-energy and the Luttinger-Ward functional. Fiy^ik = N'^'^TTt In Gq is the 
free energy of the conduction electrons. In ( ^Tj) the formulas TrlnG are ambiguous. We 
must precisely define which regularisation of these sums we consider : the actual value of 
the sums depends on the precise definition of the functional integral. For the fermionic field. 
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the standard procedure of adding and substracting the contribution of a free local fermion, 
and introducing an oscillating term to regularize the Matsubara sum holds: 

Tr In = -T In 2 + T ^ In {iuJnGf{iuJn)) e^"^""^ (48) 

n 

The situation is somewhat less familiar for the bosonic field. As pointed out above, the 
latter is merely a commuting auxiliary field (rather than a true boson). We have found that 
the correct regularisation to be used is: 

n=N 

TrlnGB = Tlim Y In (JGBiicOn)) (49) 

n=-N 

The factor of J takes into account the determinant introduced by the decoupling with B, 
and a symmetric definition of the (convergent) Matsubara sum has been used. Some details 
and justifications about these regularisations are given in Appendix (0). 
We shall perform a low-temperature expansion of Eg. (|17|) , considering successively the 
particle-hole symmetric (go = 1/2) and asymmetric (go 7^ 1/2) cases, which require rather 
different treatments. 



1. The particle-hole symmetric point qo = ^ 

In this case A = 0, so that the first term in Eq.(^7D does not contribute. Let us consider 
the last term in (^). Using the spectral representation oi Gf and the definition of Sj we 
obtain easily (for A = 0) : 

*.f*E,(.)G,(-.)=/;;d.^ (50) 



We substract the value at T = 



^(T) - vI;(T = 0) = - r du u {pfiu, T) - pf{u, T = 0)) + 2 T du (51) 

In the second term, we can replace pf by its scaling limit. So this term is of order 0(T^^^"*"^). 
We know the asymptotics of 0/ : (j)f{x) ^x^oo GiX^^f~^ + G2x'^^f~'^. (The term x^^f~'^ 
cancels due to the particle- hole symmetry). Thus, the first term in (^) is of the form : 
T2Aj+i joo^^ x(0/(x) — Gix^^f"^) (the integral is convergent). We conclude that \1/(T) = 
\I'(0) +0(T^^-f+^), so that the last term in (|T7| ) does not contribute to the zero-temperature 
entropy in the particle-hole symmetric case. 

Let us express the remaining terms in ( P7| ) as integrals over real frequencies, using the 
regularisations introduced above. As detailed in Appendix 0, this leads to the following 
expression, involving the argument of the (finite-temperature) retarded Green's functions: 



X r+oo 

fimp / du 

TT J-oo 



n,{u) (^arctan j - - j - ,ns{u) arctan 



(52) 



In these expressions np (resp. rts) are the Fermi (resp. Bose) factor. 
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At this point, it would seem that in order to perform a low-temperature expansion of the 
free-energy, one has to make a Sommerfeld expansion of the Fermi and Bose factors. This is 
not the case however, for two reasons: i) the argument of the Green's functions appearing in 
(^) are not continuous at ca; = 0, so that a linear term in T does appear (as expected from the 
non-zero value of Simp) ii) the Green's functions have an intrinsic temperature dependence, 
and the full scaling functions computed above must be used in (^21) . More precisely, when 
computing the difference fimpiT) — fimpiT = 0), the leading term is obtained by replacing 
the Green's function by their scaling form (|33|). 

These considerations lead to the following expression of the impurity entropy (per spin 
colour) Simp = Simp/N at zero-temperature for qo = \ : 

1 _ ^ _ 1 r^°° ^1 

Simp = / du (a/(u;) - af{uj = -oo)) / du a/(u;)sgn(u;) 



/ du! {ab{uj) — ab{Oj = —oo)) H — / dujab{uj)sgn{uj)—rT-^ (53) 

In this expression, a/^e denotes the arguments of the scaling functions, obtained from (p^): 



af{ijj) = arctan ~| = — arctan (^cot(7rAj) tanh — ^ (54) 
aB{Cd) = arctan = arctan (^tan(7rA/) tanh (55) 



From ([53|), we obtain with t = tauTrA 



/ 



Simp 2 /"i , /2arctant / / u\ arctan(Mt)\ arctan(ut)\ . ^. 

— — / du \ — — u arctan — H 7 ^ (56) 



7 -|- 1 TT Jo \7r(l — M^) \ \tj u J u{l — u"^) 

To perform the integration, we note that ^ (^-^ (^(1 + t'^)j^^^ = — ^^^^^2-^ and we obtain 
finally the simple expression 



f ^ 1 o 7 + 1 ['^""WT^ ln(l + M^; 
Simpiqo = - = ln2 / du (57) 

2 TT Jo [1 + U^^ 



This can also be rewritten, after a change of integration variable, as: 



SimpijiO 1/2) — T.jSimp 

iV TT 

with 



/ -2/' ^ 



1 + 7/ 12(1 + 7) 



(58) 



rx 

fix) = / Insin(M) du 
Jo 



This coincides with the large-N limit of the CFT result, Eq.(^), in the particle- hole symmetric 
case. 
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2. The general case qo ^ 1/2 



For go 7^ |) the first term in Eq.(^) also contributes to the entropy. Indeed, as shown below, 
the Lagrange multipher A(T) at the saddle-point has a term which is linear in temperature. 
This stems from a very general thermodynamic relation, which is derived by taking the 
derivative of Z with respect to go in the functional integral, leading to: 

I I \ dF 



where the average is to be understood with the action ([Tl|). At the saddle-point, we thus 
have : A = and in particular : 



dX 



T=0 



ds^ (60) 

dqo 



We shall directly use this equation in order to compute the residual entropy, by calculating 
the linear correction in T to A, and then integrating over go- This method shortcuts the 
full low-temperature expansion of the free-energy (as done in the previous section), which 
actually turns out to be quite a difficult task to perform correctly for go 7^ 1/2 ||22[. In order 
to calculate this linear correction, we shall relate A to the high-frequency behaviour of the 
fermion Green's function. As Tif{iujn) when Un — ^ ±00 we have : 

Gfizu^) = — - -A^ + o (jT^) (61) 

This shows that —A is the discontinuity of the derivative of Gfij) with respect to r at r = 0: 

/+00 
duupfiu) = -A (62) 
-00 

Let us define g^r) by: 

where a is the spectral asymmetry parameter in Eq.(^6|) (at this stage we emphasize that 
the full finite temperature, finite cutoff. Green's function G/ is considered). Equation (|62D 
can be rewritten as: 

X = aT + tanh (^|) (9.(7(0+) - dMPl) " (5r^?(0+) + OMP')) (64) 

where we have used that Gf{0~^) + Gf{(3~) = —1. Denoting by Pg{uj) the spectral function 
associated with g, we have: 
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In the scaling limit, the spectral function pg must become particle-hole symmetric (since the 
effect of the particle-hole asymmetry in this limit is entirely captured by a in (|63D ), and 
must coincide with AfT^^f~^(j)f{Lj/T; = 1/2). Hence, following the same reasoning than 
for above, the term in (^) is of order const. + 0(T'^^f~^^). Thus we have : 



df 



a 



T=0 



dA 
df 



(66) 



T=0 



where A = drg{0~^) + drg{P~) is the discontinuity of the derivative drg. A reflects the 
particle- hole asymmetry of g and thus vanishes in the scaling limit. Actually the derivative 
1^ also vanishes as T — as we now show. Consider first sending the bare cutoff to infinity 
(along with J) so as to keep the Kondo temperature fixed. In this limit A takes the form: 
A = Tf (^)- The low energy scaling limit, in which Eq. ( ^9]) holds, can be reached by 
fixing T and sending to infinity. Since g must become particle-hole symmetric in this 
limit, this implies that f{x) vanishes at small x. Hence, taking a derivative with respect to 



temperature, of A = Tf we find that |^ ^ o ~ "^^^^ finally obtain 



dsr 



imp 



dqo 



df 



—a 



(67) 



T=0 



where a{qo) is given in Eq.(^). Integrating this equation over go (taking into account as 
a boundary condition the value of Simp{qo = 1/2) obtained above), we finally derive the 
expression of the entropy: 



Simp — Simp 



1 + 7 



/ 



TT 



1 + 7. 



/ 



TT 



1 + 7 



;i-'?o) -/ 



TT 



1 + 7 



(68) 



with, as above: /(x) = /q Insin(M) du. The expression ( |68D coincides precisely with the large 
N limit of the CFT result (^ ||23|. A plot of the residual entropy and of the asymmetry 
parameter « as a function of go is displayed on Fig.^ Sj^p is maximal at go = 1/2 and 
vanishes as go ^ as expected. 

In this section, we have discovered that the spectral asymmetry parameter ("twist") a shares 
a fairly simple relation with the residual entropy, given by Eq.(|67D. These are two universal 
quantities, characteristic of the fixed-point. Remarkably, a also coincides with the term 
proportional to T in A (while A itself is non-universal, its linear term in T is). It is tempting 
to speculate that a deeper interpretation of these facts is still to be found. 



B. Internal Energy and Specific Heat 

The low-temperature behaviour of the internal energy in the large- iV limit can be obtained 
by two different methods. We shall briefly describe both since they emphasize different and 
complementary aspects of the physics. 

In the first method, we use the effective action in the form ([TT|) , before the decoupling 
with the auxiliary bosonic field Biij) is made. We thus have a quartic interaction vertex 
between the conduction electrons at the origin and the Abrikosov fermions representing the 
quasiparticles in the spin sector, which reads: 
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J Q ^ 

T7 E (//3/«- l^'^«/3)EcLci/3 (69) 
l<a,/3<Af j=l 

One can then perform a skeleton expansion of the free-energy functional in terms of the 
interacting Green's functions for the pseudo-fermions and the conduction electrons, Gjij) 
and Gc{t). The first-order (Hartree) contribution to this functional vanishes because the 
spin operator in is written in a traceless manner. The next contribution, at second 
order, yields the most singular contribution at low temperature and reads: 

AEocJ^ f dTGc{TfGf{-Tf (70) 

At the saddle-point, the interacting conduction electron Green's function is Gc{t) oc 
G f{T)GB{—T), and hence its dominant long-time behaviour is: Gc{t) ~ 1/r. Inserting 
this, together with Gf{T) ~ l/r'^^f in (|70D, we see that the leading low-temperature be- 
haviour to the energy reads AE oc ciT^^^"^^ -|- C2T^ + ■ ■ ■, and hence to the specific heat 
coefficient: 

7 > 1 : C/T-T^^z-i _ (^^-^^ 
7 = 1 : C/T-lnl/T 

7 < 1 : G/T ~ const. (71) 

which agrees with the CFT result described above. We note that there is a quite precise 
connection between this calculation and the CFT approach: the operator appearing in 
the Kondo interaction (|69|) acquires conformal dimension 2(Ac + Af) = 1 + 2Af and has 
the appropriate structure of the scalar product of a spin current with the operator Sa/s 
(transforming as the adjoint). Therefore, it is the large- version of the leading irrelevant 
operator associated with the spin sector, as described in Sec.fT|. It is satisfying that the 
leading low-T behaviour comes from the second-order contribution of this operator in this 
formalism as well. We note that for the simple Kondo model ([^), in the scaling limit, the 
analogous irrelevant operator in the flavour sector does not show up in the calculation of the 
energy in the large-N solution. We shall comment further on this point below. 
The second method to investigate the internal energy is to push the low-temperature expan- 
sion of the free-energy to higher orders. To this end, we need to compute higher-order terms 
in the expansion ( |22a| ) of the Green's functions in the scaling regime. This computation is 
detailed in Appendix |C|. 3, and leads to: 

G,{r) = A,ir^-.s, (l) + r"'9f (^) + r^'^'s?' (^) + ■ ■ • (72) 

Gs{t) = AslS-'^^-gs + ;3-'sg' (^) + p-'-'^'g'i' (^) + • ■ ■ (73) 

Let us emphasize that the exponents appearing in this expansion are not symmetric between 
the bosonic and fermionic degrees of freedom. This is because we are dealing with the Kondo 
model for which the auxiliary field (bosonic) propagator has no frequency dependence in the 
non-interacting theory. Also, the expansion given in (|72D assumes a perfectly fiat conduction 
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band in the limit of an infinite bandwith (conformal limit). Using this expansion into the 
expression (|52D of the free-energy leads to a specific heat coefficient: C/T ~ cqT^^-^"^ + 
CiT^^f~^ + C2 + ■ ■ ■• The coefficient Cq actually vanishes, so that the behavior in ( |7lD is 
recovered. The vanishing of Cq was clear in the first approach, where it followed from the 
absence of Hartree terms. In the CFT approach, it is associated with the fact that the leading 
irrelevant operator does not contribute to the free energy at first order. The vanishing of cq 
implies non trivial sum rules relating the scaling functions gf^B and (which we have not 
attempted to check explicitly). 

We also note that this behaviour of the specific heat is modified when an Anderson model 
version of the present model is considered (as in Ref. [0). Because the non-interacting 
slave boson propagator has a frequency dependence, the exponents of the second-order 
terms as written in (^) are only correct for 7 > 1 for the Anderson model. For 7 < 
1, the term (3~^^f gf\T / 13) is replaced by (3~^ gf\T / (3) ^ while I3~^ g^'si'T / 13) is replaced by 

l^-^^'gfir/P). Asaresuh, one finds a diverging specific heat coefficient in both cases, with 
C/T ~ T~(^-i)/(T+^) for 7 > 1 and C/T ~ T-^^-^^/^^+i) for 7 < 1. The behaviour for 7 < 1 
is due to the leading irrelevant operator in the flavor sector. Similarly, for 7 < 1 in the 
Anderson model, the susceptibility associated with the flavour (channel) sector Xf is found 
to diverge [0, so that a finite Wilson ratio can still be defined as Txf/C for 7 < 1. 



C. Resistivity and T-matrix 

In order to discuss transport properties, we define a scattering T-matrix for the conduction 
electrons in the usual manner (for a single impurity): 

G{k, k', u + iO+) = Go(A?, w + i0+)5^ p + G^ik, u + iO+)r(cj)Go(^', ^ + «0+) (74) 

where G and Go denote the interacting and non-interacting conduction-electron Green's 
functions, respectively. Taking a fiat particle-hole symmetric band for the conduction elec- 
tron and denoting by po the local non-interacting density of states, this yields the local 
conduction electron Green's function in the form: 

C{uj + iO+) = J2 G{k, k') = -iTrpo (1 - inpoT{uj)) (75) 

k,k' 

Following Ref. 0, we parametrize the zero- frequency limit of the T-matrix in terms of a 
scattering amplitude as: 

Tiu = 0) ^ --^ (1 - S') (76) 
SO that, the zero frequency electron Green's function reads: 

G{^0+) = -mpo^^ (77) 

S"^ = 1 corresponds to no scattering at all, while = —1 corresponds to maximal unitary 
scattering {i.e 7i/2 phase shift and vanishing conduction electron density of states at the 
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impurity site). In the overscreened noted in [^], 5*^ is in general such that \S^\ < 1, 

reflecting the non-Fermi liquid nature of the model, and the fact that the actual quasipar- 
ticles bear no resemblance to the original electrons. In addition here, we shall find a new 
feature: namely that is in fact a complex number for non particle-hole symmetric spin 
representations {i.e qo 1/2). 

We first derive an expression for for arbitrary A^, K and spin representation Q = Nqo by 
generalizing to SU{N) the CFT approach of Ref. There, it was shown that can be 
expressed as a ratio of elements of the modular S-matrix Sa^/s of the SU{N)k algebra. De- 
noting by the identity representation, by F the fundamental representation (corresponding 
to a Young tableau with a single box) and by R the representation in which the impurity 
lives (Young tableau with a single column of Q boxes), one has 

S' = (78) 

The evaluation of these elements of the modular S-matrix can be done along the same lines 
than the conformal field theory calculation of the entropy, described above. Some details 
are given in Appendix |B[ The result is: 

(^) (-'^) - (to) (-i^^™^) 

Notice that has both real and imaginary parts in the absence of particle hole symmetry 
Qo ^ i 

Let us take the large-N limit of this expression, with K/N = 7 fixed. This reads, to first 
non-trivial order: 



Tl COS 

cot 



7r(l-2go) 
1+7 



1 + 7 sin (^) 



Ar(l+7) 



2go ^ 



sm' 



sm 



(1+7) 



^0) 



We now show how to recover this expression from an analysis of the integral equations of 
the direct large-N solution. Coupling an external source to the conduction electrons in the 
functional integral formulation of the model, it is easily seen that the conduction electron 
T- matrix is given, in the large-N limit, by: 

nuj) = l^(^ + zO+) (81) 
where Q denotes the convolution of the fermion and auxiliary boson Green's function: 

g{r) = Gfir)GB{-r) (82) 



Hence, we have at T = : 



S' = l- ^^(^0+) (83) 



We first make use of the scaling limit of the two Green's functions, given by (^), and obtain 
the scaling form of Q : 
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Gir) = Gj{t)Gb{P -t)= ""^/f ^ + • • • (84) 

(3 cosh ^ sm 

We note that, in this scaling limit, the particle-hole asymmetry of the impurity Green's 
function has been lost altogether : a has cancelled completely in the r dependence of Q 
in this limit and is only present in the prefactor. Thus, only Re can be extracted from 
the scaling limit, while ImS*^ requires a more sophisticated analysis. Eq. (|8^ ) implies 
Im^ = AjAbtt, and hence ReS^ — 1 = 2TT'^AfABPo/{Ncosh'^^). We make use of the 
expression derived in Appendix y for the product of amplitudes AjAb and obtain : 

ReS^-l = —, — ^^RetanfvrA.-— ) (85) 
(1 + 7)A^ V ^ 27 ^ ' 

After expressing a in terms of go using (^6|) as : 

tanh (I) = - cot(vrA,) tan (^^^f^) (86) 
ReS"^ coincides with the real part of (pO|). 

We now consider ImS"^, for which we need to go beyond the scaling limit and use global 
properties of the Green's functions. First expressing as a convolution on the imaginary 
axis and using dujGQ{iuj) —2i7cpo6{u!) in the limit of a flat particle-hole symmetric band 
we obtain : 



tpo-f2rrg{ii^) =A + }3{iy) (87) 



with the definitions 



A = J duGf{iu)d^Ef{iu) (88) 
B{iy) = j duj [Gfiiuj + iv)d^^f{iuj) - G f{iu)d^^f{iuj)] (89) 

In the limit of vanishing v,B{v) can be calculated from the scaling limit of Gf. We obtain : 

B{v ^ 0+) = ^— { (-e "'1+;''"' + + iit\ (90) 

^ ' 1 + 7 V ^sin27rA/ / ^ ' 

On the other hand, A contains high frequency information that is lost in the scaling limit. 
We find : 

A=-^{l-2q,) (91) 
1 + 7 

The details of these calculations are provided in Appendix 0. 



Combining (|9T| , |90| , |87| , |83D we find agreement with the large N limit of 5*1 in Eq. (pOl). 
For a dilute array of impurities (of concentration Uimp), the conduction electron-self energy is 
given by Ti{uj-\-iO^) ~ nimpT{uj), to lowest order in n^mp- As shown above, T{uj) is given in the 
large-N approach by the Fourier transform of Gf{T)GB{—T). The expansion ( [72| ) yields the 
long-time behaviour ^/(r) ~ Af /t^^s +Af /r^^f + ■ ■ ■ and Gb{t) ~ Ab/t^^^ +A%^ /t+- ■ -. 
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From the fact that 2A/ + 2Aij = 1, this imphes: G/(r)GB(r) ~ l/r + l/r^+^^z H . Hence 

the resistivity behaves as: 

p(T) ~ n,„,p„ (^^1^ - cT^^/ + ■ ■ (92) 

Where pu is the impurity resistivity in the unitary hmit. For the same reasons as above, the 
Anderson model result would lead to an exponent 2As in the regime 7 < 1 0. 



D. The limit of a large number of channels (7 ^ 00) 



We finally emphasize that all the expressions derived above greatly simplify in the limit of a 
large number of channels 7 ^ cxd. This is expected, since in this limit the non- Fermi liquid 
intermediate coupling fixed point becomes perturbatively accessible from the weak-coupling 
one [H,^. The physics of the fixed point can be viewed as an almost free spin of "size" 
Q = Nqo weakly coupled to the conduction electrons. Indeed the large-7 expansion of the 
entropy (|53D , the Green function Gf (|33|) and the twist a (|46|) are : 



s.. 



imp 



- (go In go + (1 - go) ln(l - go)) - 



7r^go(l - go) 
672 



+ ■ 



9f[r) 



cosh I 



1 + - In 

7 Sm TIT 



- + 



a 



In- 



1 - go 



+ 



(93) 



and the leading terms in these expansions are given by the corresponding quantities for a 
free spin of "size" A^go- Moreover the scattering matrix and resistivity have the following 
expansion : 



1 



Re^^ 

p(T = Q)/{nimppu) 

while the anomalous dimensions read A5 



2vr2go(l - go) 



Im 51 = O 



,-yO 



1 ,7r^ 

^go(i-go):^ + 



(94) 



2A. 



7 7" 



+ . . . , 2A 



B 



7 7' 



+ . . .. 



VII. CONCLUSION 

In this paper, we have focused on the non-Fermi liquid overscreened regime of the SU(A^) x 
Sl]{K) multichannel Kondo model. This model has actually a wider range of possible 
behaviour, which become apparent when other kinds of representations of the impurity 



spin are considered. In a recent short paper [10[, two of us have studied fully symmetric 



representations corresponding to Young tableaus with a single line of P boxes. (This amounts 
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to consider Schwinger bosons in place of the Abrikosov fermions used in the present work). 
It was demonstrated that, in that case, a transition occurs as a function of the "size" P of 
the impurity spin, from overscreening (for P < K) to underscreening (for P > K), with 
an exactly screened point in between (P = K). The large-N analysis of the overscreened 
regime P < i^' is essentially identical to that presented in the present paper for antisymmetric 
representations. 

Obviously, an interesting open problem is to understand the physics of the model for more 
general impurity spin representations, involving both "bosonic" and "fermionic" degrees of 
freedom (corresponding respectively to the horizontal and vertical directions in the associ- 
ated Young tableau). CFT methods are a precious guide in achieving this goal. In particular, 
the formulas and rules given in Appendices and ^ allow for an easy derivation of the impu- 
rity T = residual entropy and zero-frequency T-matrix, using Affleck and Ludwig's fusion 
principle and the identification of these quantities in terms of modular S-matrices. 
An open question which certainly deserves further studies is to identify which of these more 
general spin representations are such that a direct large-N solution of the model can be 
found. This question has obvious potential applications to the multi-impurity problem and 
Kondo lattice models. 

During the course of this study, we learned of a work by A.Jerez, N.Andrei and G.Zarand on 
the same model using the Bethe Ansatz method. Our results and conclusions agree when a 
comparison is possible (in particular for the impurity residual entropy and low-temperature 
behaviour of physical quantities). 
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APPENDIX A: THE STRONG COUPLING STATE 

We now describe in more details the proof of the statements in Sec.(|I|) about the nature 
and degeneracy of the strong-coupling state Rgc- For a general reference on the group theory 
material used in this appendix, the reader is referred e.g to Ref. ||18|. Let us note A/y the 



number of electrons brought on the impurity site and by Y the Young tableau with A/y 
boxes associated with the representation in which the conduction electrons on the impurity 
site combine. Because of the Pauli principle, the length of any of its hnes must be smaller 
than K (and hence A/y must be smaller than NK). Indeed, we must antisymmetrise the 
wave function separately for each flavour. 
The Kondo energy is given by : 

E = JKY.SapS'p. (Al) 

a/3 
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with 



Sc.i3 = fjp-j^S^p S'^p = clcp (A2) 

in which / denotes the pseudo-fermion and c the conduction electrons at the impurity site. 
We can introduce the hnear combinations : 

= a > p = — —j^ — a< (3 T^p = S^a a = (3 (A3) 



such that: 



This leads to the following expression of the Kondo energy: 

2E 

-r = C2{Rsc) - C^iY) - C2{R) (A5) 

in which C2{Z) denotes the quadratic Casimir operator of the representation Z. The rep- 
resentation is the specific component ofY^R associated with the bound-state formed 
by the impurity spin and the conduction electrons at strong-coupling. We recall that i? is a 
column of length Q in this paper. 

We have to minimize E over all possible choices of Y and of Rgc- 
First let us recall that for a general representation Y, C2 is given by : 

1 Tl 

C2 = -'nA-\- + l) (A6) 
where ?7,j(l<2<A^ — l)is the number of columns with length i in the Young tableau Y 



and A is the Cartan matrix of the SU(N) group [T^. Let us denote by fj (1 < j < A^) the 
length of the line j in the tableau. Then we have 

with A/y = Y,f=i fj is the number of boxes of Y. Note that with this definition, all /,'s 
can be shifted by the same constant without changing the representation (this is because a 
column of length A^ can be removed without changing the representation). ( [A^ ) can be given 
a simple interpretation in terms of A^ "particles" occupying a set of fermionic levels. This 
interpretation was introduced in a slightly different form by Douglas |12|. Let pj = fj —j + N 
be the position of the particle j. Because F is a Young tableau, the particles are ordered and 
cannot be on the same level. Fig.(|^) gives an example of the construction of the diagram 
associated with a simple Young tableau. 

A simple construction of all allowed Young tableaus appearing in the tensor product Y ® R 
|T8| can be given in this fermionic language. Starting with the diagram associated with 
Y, we choose Q particles and raise each of them by one level beginning with the one in 
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the highest level. (We note that, in the fermionic interpretation, adding a box to line i 
corresponds to raising the i-th particle up by one level). An example is given in Fig.(|]). 
Let us denote by Pi the positions of the "particles" in Y, and by the new positions in 
a given allowed component of F ® i?. The Kondo energy is given by: 

^ = ^ I UpI - Pi) - § - ^ - c.{R) I (A8) 

The last two terms are constant {R is held fixed) and can be dropped in the minimisation 
process. 

The Pi's can be decomposed in two sets: those for which p[ = pi (we have N — Q oi them) 
and those for which p'i = Pi + ^ {Q of them). Let us denote by P the sum of the latter ones. 
We have: 

Thus the lowest energy is achieved for the smallest possible value of P. Since a given shift 
p ^ p + 1 can only appear once in the sum (because double occupancies are forbidden and 
a given particle cannot be raised twice), the absolute minimum is obtained when we sum on 
all the lowest Q shifts. This implies that the diagram associated with Y has Q particles on 
the Q lowest levels (from to Q — 1) and none on the Q-th level. 

The upper part of the diagram (above level Q) is then determined by the maximisation of 
A/y. Going back to the language of Young tab leans, the minimum is thus achieved when Y 
is a rectangle of height N — Q and width K, and Rgc is given by the same tableau with the 
first column removed. 
Two cases must thus be distinguished: 

• For [K = 1) and for arbitrary A^ and Q, Rsc is the trivial (singlet) representation (of 
dimension d{Rsc) = !)• 

• For [K > 2, arbitrary A^ and Q) the dimension d{Rsc) is is larger than the dimension 
of R. Indeed, denoting by dxiRsc) the dimension Rsc for K-channels, we have the 
recursion relation (from the "hook law" pO[]): 



dK+i ^ {N + K){N + K-1)...{N + K-Q + 1) ^ 
dK {Q + K){Q + K -1)---{Q + K + 1-Q) ^ ' 

(because Q < N). It increases with K. The K = 2 case is just a column of length 
N — Q which has the dimension of R. Moreover the inequality is strict for K > 2. 



APPENDIX B: CONSTRUCTION OF MODULAR S-MATRICES 

If two representations R and R' correspond to fermion configurations with positions 
{pi, . . . ,Pn} and {p[, . . . ,p'j^} (See Appendix |^) respectively, then the modular ^-matrix 
element is : 
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s 



R,R' 



27r iNpp 



det \e N+K 



(Bl) 



with p = J2iPi/N, P' = J2jP'j/N and C^^k is a constant which depends only on and K. 
Since for the trivial representation 0, the p's are the consequent integers 0, 1, ... , N—2, N—1, 
So,R involves a determinant of the form 



1 


1 • 


•• 1 


1 




Z2 ■ 


■■ ZN-1 


zn 




Z2 


■ ■ ^2 

ZN-1 


Zn 


N-2 
1 

Af-1 
1 


^N-2 
Z2 
N-1 . 


^N-2 

zn-1 


^N- 
Zn 

Zn 



2 Trip J- 

where Zj = , {pi, . . . ,Pn} being the positions of fermions corresponding to the repre- 
sentation R. This is just the Van Der Monde determinant A{z) = Yli^j{zi — zj) 
To calculate the T matrix, we also need to know the ^-matrix element between the fun- 
damental representation F and an arbitrary representation R. For F, the positions of the 
fermions are 0, 1, 2, — 2, N. Therefore Spu involves the determinant 



1 



^N-2 



1 

Z2 



1 



1 



zn-1 Zn 



z 



Z2 



N-1 



^N-2 
^N-1 
^N 
^N-1 



Z 



N 



,N- 
-N 



,N 
-N 



This determinant has the same antisymmetry property in ZiS as the Van Der Monde deter- 
minant. However, the present determinant is one order higher than A (2;) as a homogeneous 
polynomial in z's. A little reflection shows it is {J2^Zi)^{z). 
Finally we find : 



S 



F,R 



So,. 



27vip 



N 



(B2) 



Using this formula we deduce (|T^) from (|78| 



APPENDIX C: SOLUTION OF THE SADDLE POINT EQUATIONS IN THE 

SCALING REGIME 

In this appendix we solve (|1^) in the scaling regime as explained in section (|V A| ) and obtain 
the scaled spectral densities and Green functions. 
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1. Scaling functions 

First we show that ( PPD is solution of the saddle point equations in the scahng regime. We 
deal with an arbitrary go- Let us denote by ct/^b the scaling function of the fermionic and 
bosonic self-energies : 

S/,B(r)=A^,//?~'^-'^-'V;,B(^) (CI) 

Go is the local Green function for the conduction electron. Its density of states does not 
depend on T. So its scaling form is : 

with pq = — i ImGo(i^ = 0). Using this formula and (|15|), ct/ b are related to gf^B- We insert 
the scaling form ( p2a| ) into (|16D. Matching the power in (3 leads to 2Aj + 2Ab = 1 and : 

Aj'^gJ^{iUn) = (A - T^fiiuo)) P^^^ - Ab {(yf{iu)n) - crf{i^^o)) , 

AB^QB^^^n) = {j- Sb(0)) P'^'f - Af {aB{^Un) - (Tsm (C3) 

with = iC^n + l)7r and zz/„ = i2n7c. The term iUn in (|T^) vanishes in this scaling limit 
because Af^B < 1- We assume that at zero-temperature 

A - S;(0) = ^ - S^(0) = (C4) 
so P disappears of these equations at lower-order. 

We insert our Ansatz into ( |C3|) with the following Fourier transform formulas (which follow 



from 111 3.631) 



(27r)^^/T^-^/-^.(-l)»+^r(l - 2A;) 

qriiUn) = 7 1 X 7 — r Coa 

coshfr(l-A,-|^ + f)r(l-A, + |^-^) 



(2vr)^^^T^^^-H-l)"+^r(l - 2A^) 
gB{tJ^n) = 7 1 — 7 : (C5b) 

coshfr(i-A^-|^ + f)r(i-A^ + |j 



27r 27r 



afilOJn) — 0-f[lUJn) = 7 r 7 r (u = 0) C5c 

' ' coshfr(i-A,-i^ + f)r(i-A, + |^-f) 



-shfr(i-A,-|j + f)r(i-A, + |^-f) 



We see that (pOD is solution of (|C^ ) provided that : 
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The precise form of the cancelation 
T): 



-^-Sb(0) = - , 

Jk cosh f r (i 



at finite temperature is (at leading order in 
PoA/(27r)2^/+iT2^/r(-2A/) 



A 



A-S 



z7PoAB(27r)^^^+iT^^^r(-2AB) 

coshfr (-A^ + f)r(i-A^-f) 



(C6) 
(C7) 



(0) is obeyed : 2A; = ^ , 2Ab = ^ 
We have the relation between amplitudes : 

2Ab = -2-iAfABPoT{l - 2^b)T{2^b) 



sin (^ttAb — Y 



cosh 



(C8) 



In our scaling forms, a is the same for the fermionic and the bosonic function. One can 
check easily that for a more general Ansatz with af and as the saddle-point equations imply 
af = Ob- 



2. Spectral densities 

We calculate now the scaled spectral density from the above scaling function. Denoting 
C = — 1 for fermions and C = 1 fo^' bosons we have the general formula. 

_ l_t^e-Ps PfM^)de 0<r<(3 (C9) 

In the scaling regime, we have to solve : 



/ + 00 g — J-Ii 

_J f,B{^)du 0<x<l (CIO) 
-oo i — C 6 



cosh I \sm{nx) J J-oo 1 — 

Setting t = i(x — |) we see it is sufficient to solve : 



/ — -—u(f)f B{u)du |Imt|<- (Cll) 

j-oo e 2 — Ce 2 2 



cosh I ycOsh(7rt)y J-oo 62—^6 

Due to the properties of Fourier transformation, we can just solve for a = 0, and obtain the 



solution for arbitrary a with (0). With 

, A 



TT 



(jji " \ ^-itu _ {2 + Z)^ {2 Z) ro<A<i 

00 \cosh(7rt) I r(A) I u real 



(see formula 3.313.2 of [^]), we find the result given in the text 
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M^, % = 1/2) = i(2vr)2^/-i cosh f "^^^^V|a,)^"^^^ 

M^, % = 1/2) = ^(27r)2^-i sinh f "^^"%f I^'^""^^^ (C13) 

The asymptotic behaviour follows from formula 8.328 of 

We then derive the full Green function by taking a Hilbert transform : 




(CM) 



We find ( |33D using : 
• The representation 

= r°°e^^("-")dA lmz>0 (C15) 
z — u Jo 



The Fourier formula which inverses (|C12| ) 
The formula 



dx- 



ifSz 
2tt 



r(i-A) 



(^sinh ^ 



TT 



< A < 1 
2; real 



(C16) 



which results from formula 3.112.1 of ||21||. 



Finally, we comment on the treatment of the constraint equation (j[7\) in our derivation of the 
scaling funtions. The relation between a and go has been derived from a Luttinger sum-rule, 
which holds at zero-temperature. So one may worry whether the scaling form does satisfy 
the leading low-temperature corrections to the T = constraint equation. We show now 
that this is actually the case. Starting from Eq.(p!7D written as: 



(C17) 



we substract the relation at T 
given by Eq. (pSj) to obtain : 



dx 4>f{x) — 



and take into account the asymptotic behaviour of 



cosh fF (2 A 




°° , sgnxchfix) 
dx — -—^ 

el^l + 1 







(C18) 



It is a rather strong constraint on the scaling function 
and it is satisfying that the explicit form obtained for 
that 0/ is really a solution of the full system (15,16,17) in the scaling regime at fixed go 



that this equation should hold, 
does satisfy ( |C18D . This proves 
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3. Higher-order terms in the scahng expansion 

Here, we give some indications on the derivation the expansion in (^). Let us start from 
the long-time expansion for T^^ <^t <^ (3: 

G/W~;^ + ^ + ^ (C19) 

in which a and A are exponents to be determined below. Then we have : 

Gf{uj) ~ A^CaA,-!^'""^'"' + (C20) 

it 

with Ca = / rft^s+T, and a similar expression for Gb- We can then deduce the expansions 
of and S^, and insert them into the saddle point equation. We find : 



■uo • ^ = U) U) ° UJ 



The first order yields: 

-71 = PoAfABC2AfC2AB-l = 7P0^/^bC'2AsC'2A/-1 (C22) 

which, using: 

Ca-1 oc ACa (C23) 



gives (|19D again. The second equation leads to A = a + 1 — 4Aj. 

First suppose A < 1 : in this case we must drop the u term but we have 

Ca-lC\-i CaC\ 



C2AB-1C2AJ-I C2AbC2A 



(C24) 



which implies a = 2Af or a = 2Af — 1 (taking |C23| into account). So this possibility must 
be rejected. Finally we are lead to A = 1 and a = 4Ay. 

The higher order corrections can be dealt with in a similar manner. Restoring the scaling 



functions, this leads to ([721). 



APPENDIX D: CALCULATION OF THE RESIDUAL ENTROPY 

1. The formula of the free energy 

We first give a few more details on the regular isat ion in (^7]). We will check that (^) is the 
right formula for the pseudo-boson. 

In the following we will denote by Tr-t the regularisation with 6*"^"° and by Tigym the 
regularisation defined in (|49|). We note that Tr^^^ = (Tr^ + Tr_)/2 as can be checked 
explicitly using a spectral representation of the function to be summed. 
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Let us introduce the following notation for any quantity A (function of A) : A\A = A^ — A ^. 
As the free energy is particle-hole symmetric, we have : 



A = AA(rTrlnG/) -7AA(rTrlnGB) (Dl) 



Let us consider 



ILJ„ + A — Tij^iuJn) 



such as 



(toOn) = In — ^ (D2) 



AA(rTr+lnG'/) = -</)(r = 0") (D3) 



As (f) is particle-hole symmetric, we have (j){T = 0+) = — 0(r = 0^). As its asymptotic 
behaviour is (j){iiOn) ~ its discontinuity is 0(r = 0^) — 0(r = 0^) = — 2A. 
We obtain : 

AA(TTr+lnG/) = -A (D4) 



This implies that the bosonic term does not contribute to (pi|) . But there is an analogous 
relation for the boson : we first calculate the discontinuity of from the saddle-point 
equations, use an analogous function (p and obtain AA(TTr-|- In Gb) = t(1 — 2go)^/2. So we 
find 

AA(TTr,j,^lnGB) =0 (D5) 
So we have checked that (HO) is the right regularisation for the bosonic term. 



2. Derivation of ( |5^) 

We consider first the fermionic term. Let GniiUr,) = ^ be the Green function of free 



electrons. We have 



TTr+lnGf = -Tln2 - - / dcu (ImlnG f - ImlnGo) npiuj) 

■K Jn 

= -T\n2 + - [ duj f arctan§^4 + ^ - ^^(-^) I ^i^M (D6) 

'K Jn \ G f[x) 2 J 

1 / G'Ax) 7r\ 

= - / duj\ arctan f - - npiuj) (D7) 

7T J— oo y (j j\X) Z J 



The bosonic term is obtained by an analogous calculation. In the particle hole symmetric 
case considered in the text the three regularisations for the bosonic term are equivalent. We 
have : 



\ r+oo 

— TTr^ym IiiGb = / dLjlm\n{JG{Lj))nB{uj) 

1 /•+^ , G'Ux) , , , , 

= / acjarctan nB\uj) (D8) 

TT J-oo G nix) 



Finally we find the formula quoted in the text Eq. (^). 
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APPENDIX E: SOME DETAILS OF THE T-MATRIX CALCULATION 

In this appendix, we calculate A and B{0^). 

1. Computation of A 

Using the definition of S and introducing a oscillating term to regulate the 2 integrals, we 
have : 

A = J dujGf{iu)d^'Lf{iu) (El) 
= i r dijGf{iu})e'^^^ + r dujdJnGf{iuj)e'''^^ (E2) 

J —OO J —QO 

The first term is 2mqQ. Using (^) and 



duj = ZTT (E3) 

OO UO 

(the integral is to be understood as a principal part), we have with iIj{z) = ln{zGf{z)) 

A = -2ie{l + 7) + / dz{d,^){z) (E4) 

JiTZ 

= -2ie{l + 7) - lim 2i lm{ip{ie) - ip{ioo)) (E5) 
= -2^7^ (E6) 

We used ip{ioo) = and ip{ie) ~ \ia{Ae^^f) — 16 with A a real constant. Finally we find (^TJ). 



2. Computation of ^(0+) 
B{v) = i I du {Gfiiu + iv) - Gfiiu))) - [ duj {Gfiiio + iv) - Gfiiu)) d^Gfiiu) (E7) 



We replace G/ by the scaling function gj. The second term is of order 1 whereas the first is 
OiT'^'^f) and can be neglected. We have then 



= - J dx {gf{i{x + l)z/) - gjiixv)) d^g/iixu) (E8) 

with 9 = j;. We want i3(z/ = 0~^,T = 0) which is obtained by taking the limit z> +00 in 
the previous scaling limit of B. To perform this limit we use gf{z) = gf{z) and the following 
expansion for g : 

g f{ix) ^x-*+oo cAx'^^f^^ with A = icosh(^— + inAj^ (E9) 
where c is a real constant (Eq. (E9|) is obtained directly from (p3D). We find : 
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S(0+) = (1 - 2A;) 



dx- 



-1 

oo 
ixO+ 



dx- 



-=■ [ dX; 



— / dx 



+ 



X 



(ElO) 



The last term a principal part and is given by ([E3|) . We then use the following identity : 







-1 

-oo 
oo 



dx- 



+ e-^'^^f I dx 



/•oo 

+ / dx 
Jo 



rO+ 



1 |a;|2^/(x + 1)1-2^/ 



2%(a; + 1)1-2^/ 



We find 



i5(0+) = (1 - 2A^) 



A 

3 + ^ 



~2i7rA 



/'^|X|2A/(X + 1)1-2^/ 



7 



1 + 7 



TT 



sin 27rA 



«7r 



/ 



A simple calculation with (j^^ shows : 

A _ 
A ~ ' 



.g-T^(l-2go) 



(Ell) 



(E12) 



(E13) 



and we find flOOD. 
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FIG. 1. Young tableau corresponding to the strong-coupling state 
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FIG. 2. Plot of ^/ and (f)^ as a function of a)— a for different values of the asymmetry parameter 
a : a = 0, —2, —5, — oo(Ay = .3) 



37 




38 




FIG. 5. An example of the general composition rule explained in the text. We show the 
fermionic diagram associated with Y, the resulting fermionic diagrams and their transcription in 
terms of Young tableaus. {N is arbitrary in this example as evidenced by the ) 
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